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Abstract. Efforts to understand the deviation of the L-T relation from 
a simple scaling law valid for clusters and groups have triggered a num- 
ber of interesting studies on the subject. Techniques and approaches 
differ widely but most works agree on the important role played by gas 
cooling and heating sources like AGNs and SNe. Observations set use- 
ful constraints on the evolution of the intracluster medium (ICM): a 
lOOKeV/cm^ entropy floor in the core of groups and about 5-15% of 
baryons being converted into stars. However, essential details like the 
nature of the dominant heating mechanism and the quantitative impor- 
tance of cooling still need to be addressed. I suggest that a new generation 
of high resolution N-body simulations and a quantitative comparison of 
results between different approaches is required to improve results and 
increase our understanding of the problem. 



1. Q: What are the Theoretical Expectations for the L T relation? 

The first attempt to model the ICM in the framework of the hierarchical scenario 
assumed its thermodynamical properties to be entirely determined by gravita- 
tional processes, like adiabatic compression during collapse and shock heating 
(Kaiser 1986). If no characteristic scales are present in the underlying cosmology 
(i.e., Einstein-de-Sitter cosmology and power-law shape for the power spectrum 
of density perturbations), this model should predict hot gas within rich clus- 
ters to look the same as within poor groups, since gravity in itself does not 
have characteristic scales. Under the assumptions of emissivity dominated by 
free-free bremsstrahlung and of hydrostatic equilibrium of the gas, this model 
predicts Lx oc r^(l -|- z)^/^ for the shape and evolution of the relation between 
X-ray luminosity and ICM temperature (also Eke, Navarro & Frenk 1998). 

Furthermore, if we define the gas entropy as = T/ue^^ {Uf,: electron number 
density; e.g., Eke et al. 1998), then the self-similar ICM has S oc T(l -|- z). 

2. Q: What are the observational Constraints? 

However, soon the simple model described above failed to account for several 
new observational facts: (a) the Lx-T relation for nearby clusters is steeper 
than predicted, with Lx oc for T > 2 keV clusters (e.g., David et al. 1993; 
White, Jones & Forman 1997; Allen & Fabian 1998; Markevitch 1998); with 
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a possible further steepening at the group scale, T < 1 keV (Penman et al. 
1996; Helsdon &: Penman 2000); (b) no evidence of evolution for its amplitude 
has been detected out to z > 1 (e.g., Fair ley et al. 2000; Delia Ceca et al. 
2001; Borgani et al. 2001a); (c) the gas density profiles in central regions of 
cooler groups is relatively softer than for hotter cluster and, correspondingly, 
the entropy is higher (e.g., Ponman, Cannon & Navarro 1999; Lloyd-Davis et 
al. 2000). Observational evidence provides a fairly robust measure of the entropy 
floor for the diffuse hot gas in groups (Lloyd-Davis et al. 2000) of 100 KeV/cm'^. 
A related, very important measurement is the fraction of baryons that has been 
converted into a "cold", non X-ray emitting phase (i.e stars and neutral hy- 
drogen). This constraint is crucial as it defines the amount of energy available 
from SNe and the importance of cooling in groups and clusters (Renzini 2000). 
Balogh et al. (2001) estimate that only a small fraction (5-10%) of baryons has 
been converted into stars, almost independently of the total virial mass of the 
cluster. 

3. Q: What are the physical processes that could affect the L— T 
relation? 

• Gas cooling. At the center of groups the cooling time is much shorter than 

a Hubble time. Cooling would remove low entropy gas transforming it into 
stars and possibly originate the entropy floor observed. While advocated on the 
ground of its simplicity (Muanwong et al 2001 also Voit & Bryan 2001), the 
amount of gas involved cannot exceed the observed fraction of 5-15% of the 
total baryon fraction. 

• Energy injection from AGNs (e.g Bohringer, this proceedings) would be an 
important source of heating. Most of the energy would come from jets that 
would carve "bubbles" in the ICM (Mc Namara et al 2001, Quilis et al. 2001), 
while the X-ray emission from the central engine in radio quiet QSOs would 
likely have a small effect due to the small cross section for the process. 

• Energy injection from, SNe. Winds propelled by SNc explosions (an hefty 
lO'^'^crg/SN) would transfer heat to the ICM, possibly escaping the halos of 
individual galaxies in the forming protoclustcr and enrich the ICM with metals. 
(Note: substantial metal enrichment could originate also from ram pressure and 
tidal stripping of galaxies (Moore, Quilis & Bower 2001). 

4. Q: What are the most used theoretical tools to study the L— T 
relation? 

• Analytical methods allow a fast exploration of parameter space, but must some- 
times rely on the parametrization of the end results of complicated physical pro- 
cesses. While this approach requires studying the physics of gas in a somewhat 
idealized situation (spherical collapse, accretion from a uniform background) 
the amount of insight gained is impressive, although conclusions sometimes dif- 
fer widely (Kaiser 1991,Valageas & Silk 1999, Tozzi & Norman 2000, Bryan 
2000, Balogh, Babul & Patton 1999, among many). The most refined methods 
are based on the standard machinery developed to link gravitational clustering 
in Cold Dark Matter Models and the physics of gas and star formation (Cole 
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et al. 2000, Somerville & Primack 1999, Bower et al. 2001, Menci & Cavaliere 

2000, Cavaliere, Giacconi & Menci 2000). 

•N-body simulations are a powerful method, able to treat the highly non-linear 
problem of the formation of cosmic structures with very few assumptions. Recent 
observational evidence for a mono phase ICM (Bohringer et al 2001) remove 
some potential worries about the use of SPH methods, which do not describe 
well a multiphase medium. Some recent works on the L-T relation including 
feedback and star formation (Mctzler &: Evrard 1994, Bialek, Evrard &: Mohr 

2001, Borgani et al 2001b, Valdarnini 2001) provide a good starting point and 
an extensive list of references to the subject of simulations of galaxy clusters. 

5. Q: Heating, Cooling or both? 

As mentioned above, several papers suggested that simply the addition of cooling 
would remove enough gas from the ICM to reconcile models with the observed 
L-T relation. This would likely be in contrast with observational evidence that 
points to a limited fraction of baryons in clusters in a "cold" non X-ray emitting 
phase. Moreover it would be difficult to hide the large quantity of cooled gas in 
"dark baryons" as they would cause detectable deepening of the central potential 
well at the center of clusters (e.g Lewis et al. 1999, Valdarnini 2001). This 
would have a substantial effect on the predicted M-T relation increasing the 
emission weighted temperature of groups (Finoguenov et al 2001). On the other 
hand recent theoretical works, both analytical and numerical seem to (perhaps 
slowly) converge to a required energy input comparable to 100% of what is 
available from SNe (Kravtsov & Yepes 2000, Bower et al 2001, or even a few 
times higher (Borgani et al. 2001b, Valageas &i Silk 1999, Wu, Fabian Sz Nulsen 
2000) to satisfy the L-T relation. While semi-analytical models estimate the 
energy available from SNe to be of the order of ~ 0.25 keV per particle , Renzini 
(2000) suggested an even smaller energy budget available from SNe: just 0.1 keV. 
Taken together these works suggest that the required energy must be provided 
by both SNe and AGNs, perhaps active as major mergers shape the early type 
population of cluster galaxies. Unfortunately none of the N-body simulations 
cited above (some analytical works did) treated both cooling and heating in a 
fully consistent way or included the effect of AGNs on the ICM. This leaves the 
quantitative contribution of cooling and heating rather uncertain. 

6. Q: Is there a preferred epoch for heating the ICM? 

Not really. As the entropy S depends on a inverse power of p it might seem more 
efficient to heat the ICM at lower z (say z< 5) as suggested in Tozzi & Norman 
(1999). However, in their model the authors heated the gas before accretion, 
while at the background density. In a more realistic case gas will be heated 
when already inside halos that will later merge to form the cluster. While this 
requires more energy for a given entropy S, (p is higher by a factor ~ 200) this 
is more than compensated by the fact that halos collapsed at high z will form 
the core of the cluster at lower z (Governato et al. 2001), i.e. you dump energy 
right where it is needed. Borgani et al (2001b, 2001c) and Bialek, Evrard & 
Mohr (2000) tried a large range of z (1 to 5 and 20 respectively) finding a weak 
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dependence on results for the argument given above. Of course, a redshift too 
close to z ~1 for heating would cause the L-T to evolve significantly at moderate 
redshift, contrary to observations. 

Bottom line: it is very likely that a combination of cooling &: feedback will give 
the correct, non-evolving L-T relation, although only stronger observational 
constraints and more refined simulations will allow us to evaluate the quantitative 
role of cooling, SNe and AGNs. 

7. Future Improvements: 

Numerical simulations including a treatment of hydrodynamical processes have 
reached a high degree of maturity, with widely different codes giving comparable 
results (Frenk et al 1999). However, the number of particles used to simulate 
individual halos is still systematically lower than in DM only runs (e.g Ghigna 
et al. 2001). The hot debate on the shape of inner density profiles in DM 
halos has started a race to higher resolution runs (Governato, Ghigna & Moore 
2001). While the quest for the ultimate density profile of DM halos might 
ultimately prove futile it has sparked a number of improvements on codes and, 
equally important, increased the robustness of results. A key method is the so 
called "renormalization technique" (Katz Sz White 1993) where a halo taken 
from a cosmological simulation is resimulated at much higher resolution, while 
sparsely resampling the surrounding region. Fig.l shows the crucial effect of 
increasing resolution on a sample of SPH simulations of renormalized halos. 
(Borgani et al 2001b, c). The L-T relation obtained from simulations of large 
cosmological volumes (where particle number and the ratio softening/R^j^ for 
each halo depends on its virial mass) likely suffer of systematic biases due to 
the worsening resolution for less massive halos: namely smaller hx and slower 
cooling due to lower central densities. 

This is my very personal wish list for future theoretical and observational studies 

on the subject: 

• Estimates of Lx and T in groups that exclude contamination from AGNs. 

• Estimates of the amount of metals in clusters to a significant fraction of K^ir as 
another method to infer the energy release from SNe (Renzini 2000, De Grandi 
& Molendi 2001). 

• A robust estimate of the fraction of "cold" baryons which keeps into account 
the amount of intracluster stars. Recent claims (sec Arnaboldi et al 2001) based 
on intracluster PNs almost double the the efficiency of star formation in clusters. 

• Simulations of sufficient resolution to treat cooling and star formation processes 
in a robust way. (note: robust does not mean realistic, but it's a start). 

My checklist for simulating an individual cluster: 

N>10^ 

< 1% 

^gas{r <Rvir) > 10^ 

Simulation box lOOMpc 

N steps » 10^-4 
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Figure 1. X-ray luminosity vs mass weighted Temperature: effects of 
resolution. Empty (red) dots are runs with ~ 0.5-2 x 10^ gas particles 
within Hyir. Full dots are halos simulated with Ngas{r >Ilyir) ~ 10^. 
Softening is l%Rmr- The upper filled (red) dots show the effect of 
having instead e = 0.5% Rvir- Lz grows by 15-20%. (Green) filled 
squares show the large effect of including line emission in the estimate 
of X-ray luminosities for the high res runs. The dashed line is the 
scaling predicted by Kaiser (1986) assuming NFW (Navarro, Frenk & 
White 1995) profiles. ( data from Borgani et al. 2001 in prep.)Bottom 
line: use at least 1(P gas particles per halo to avoid nasty systematic 
effects like reproducing a bend in the L-T relation at T < 0.3 KeV. 
Correct for line emission if you compare vs. real data. 
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• Exploring the role of QSO's and AGNs in heating the ICM in a cosmological 
context. 

• Introducing heating schemes that closely follow the time evolution of the phys- 
ical process behind it. 

• Resolving the discrepancy between adiabatic runs with mesh and SPH codes 
(e.g Kravtsov et al. 2001, also Springel & Hernquist 2001 and Cen & Ostriker 
1999) : why SPH simulations do not have an entropy floor and mesh codes do? 

Acknowledgments. I wish to thank Stefano Borgani, Marino Mezzetti and 
Riccardo Valdarnini for organizing such an interesting and lively meeting and 
Uros Seljak and Bepi Tormen for stimulating discussions (although they made 
me skip a coffee break or two). I also thank my collaborators Stefano Borgani, 
Nicola Menci, Tom Quinn, Joachim Stadel, Paolo Tozzi & James Wadsley for 
allowing me to discuss our results and show Fig.l. in advance of publication. 

8. References 

Allen S.W., Fabian A.C., 1998, MNRAS, 297, 63 
Arnaboldi et al, 2001 AJ in press 

Balogh M.L., Babul A., Patton D.R., 1999, MNRAS, 307, 463 

Balogh et al., 2001, MNRAS, 326, 1228 

Bialek, J.J., Evrard, A.E., Mohr, J.J, 2001, ApJ, 555, 597 

Bohringer, H. et al 2001, A&A, in press 

Borgani et al. 2001a ApJ, 561, 13 

Borgani et al. 2001b, ApJ, 559, 71 

Borgani et al. 2001c, in preparation 

Bower R.G., Benson A.J., Bough C.L., Cole S., Frenk C.S., 2001 MNRAS 325, 
497 

Bryan G.L. 2000, ApJ, 544, 1 

Cavaliere A., Giacconi R., Menci N., 2000, ApJ, 528, L77 
Cen R., & Ostriker J.P., 1999 ApJ517, 31 

Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, MNRAS, 319, 168 
David L.P., Slyz A., Jones C, Forman W., Vrtilek S.D., Arnaud K.A., 1993, 
ApJ, 412, 479 

De Grandi S. & Molendi S., 2001 ApJ, 551, 153 

Delia Ceca R. et al. 2000, A&A, 353, 498 

Eke, v., Navarro J.F. & Frenk C.S. 199,8 MNRAS, 503, 569 

Fairley B.W. et al., 2000, MNRAS, 315, 669 

Finoguenov A., Reiprich T.H., Bohringer H., 2001, A&A 

Frenk, C.S. et al. 1999, ApJ, 525, 554 

Ghigna S. et al. 2000, ApJ, 544, 616 

Governato, F., Ghigna, S. & Moore, B. 2001, proceedings of "Astrophysical Ages 
and Timescales, ed. Ted Von Hippel, (|astro-ph /0105443 ) 

Governato, F., Ghigna, S., Moore, B., Quinn, T., Stadel, J., & Lake, G. 2001, 
ApJ, 547, 555 

Helsdon S.F., Ponman T.J., 2000, MNRAS, 315, 356 

Kaiser N., 1986, MNRAS, 222, 323 

Kaiser N., 1991, ApJ, 383, 104 

Katz, N. & White, S.D.M., 1993 ApJ412, 455 



APS Conf. Ser. Style 



7 



Kravtsov A.V., Yepes G., 2000, MNRA S, 318, 227 

Kravtsov A.V. et al. ApJ submitted, ( |astro-ph/0109077| ) 

Lewis G.F., Babul A., Katz N. Quinn T., Hernquist L., 2000, ApJ, 536, 623 

Lloyd-Davies E.J., Ponman T.J., Cannon D.B., 2000, MNRAS, 315, 689 

Markevitch M., 1998, ApJ, 504, 27 

Mc Namara, B.R. et al 2001, ApJ, 562, 149 

Muanwong et al 2001, ApJ2001 552, 27 

Menci N., Cavaliere A., 2000, MNRAS, 311, 50 

Metzler C.A., Evrard, A.E. 1994 ApJ437, 564 

Moore B., Quilis V. & Bower, R. 2000, ASP Conf. Sr. 197, 363 

Ponman T.J., Bourner P.D.J. , Ebeling H., Bohringer H., 1996, 293, 690 

Ponman T.J., Cannon D.B., Navarro J.F., 1999, Nature, 397, 135 



Quilis, v.. Bowe r R.G. fc Balogh, M.L. 2001, |astro-ph/0109022 

Renzini A., 2000 astro-ph /OOO 13 1^ , proceedings of "LSS in the X-ray Universe, 

1999, ed Atlantisscience. 

Somerville, R.S., & Primack, J.R. 1999, MNRAS, 3 10, 1087 
Springel, V. & Hernquist L. 2001 |astro-ph/01110T6 



Tozzi P., Norman C. 2001, ApJ, 546, 63 

Valageas P., Silk J., 1999, A&A, 3 50, 725 

Valdarnini R., 2001, ApJ, in press |astro-ph/0110545| 



Voit G.M., & Bryan G.L., 2001, Nature, in press 

White D.A., Jones C, Forman W., 1997, MNRAS, 292, 419 

Wu K.K.S., Fabian A.C., Nulsen P.E.J. , 2000, MNRAS, 318, 889 



